Disorder and interaction effects in two dimensional graphene sheets 
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The interplay between different types ol disorder and electron-electron interactions in graphene 
planes is studied by means of Renormalization Group techniques. The low temperature properties 
of the system are determined by fixed points where the strength of the interactions remains finite, 
as in one dimensional Luttinger liquids. These fixed points can be either stable (attractive), when 
the disorder is associated to topological defects in the lattice or to a random mass term, or unstable 
(repulsive) when the disorder is induced by impurities outside the graphene planes. In addition, we 
analyze mid-gap states which can arise near interfaces or vacancies. 



Introduction. Graphite is a widely studied material, 
which has attracted recent interest due to the observation 
of anomalous properties, such as magnetism or insulating 
behavior in the directionperpendicular to the planes in 
different samples |1I1I1I1I1IEI3,I1. 

The conduction band of graphite is well described by 
tight binding models which include only the tt orbitals 
which are perpendicular to the graphite planes at each 
C atom0. If the interplane hopping is neglected, this 
model describes a semi metal, with zero density of states 
at the Fermi energy, and where the Fermi surface is re- 
duced to two inequivalent points in the Brillouin Zone. 
The states near these Fermi points can be described by 
a continuum model which reduces to the Dirac equation 
in two dimensions. Due the the vanishing of the den- 
sity of states at the Fermi level, the long range Coulomb 
interaction is imperfectly screened. This implies that a 
standard perturbative treatment leads to logarithmic di- 
vergences, and to non trivial deviations from Fermi liq- 
uid theory[l3 IH 113 ■ 

In the strong coupling regime, 
the model can exhibit a phase transition which leads to 
a rearrangement of the charges and spins within the unit 
cell, and which is similar to the chiral symmetry breaking 
transition found in field theories |13lll4| . 

It is known that disorder significantly changes 
the states described by the two dimensional Dirac 
equation[ll [H [l3|, and, usually, the density of states 
at low energies is increased. Lattice defects, such as pen- 
tagons and heptagons, or dislocations, can be included 
by means of a non Abelian gauge field 0, 0|. In gen- 
eral, disorder enhances the effect of the interactions. In 
addition, a grap hene plane can show states localized at 
interfaces pOl l2ll | . which, in the absence of other types of 
disorder, lie at the Fermi energy. Changes in the local 
coordination can also lead to localized states 22). 

The present work attempts to study, within the same 
footing, the role of long range interactions and disorder. 
This problem has already been studied in relation with 
critical points between integer and fractional fillings in 
the Quantum Hall Effect [23L |23 |. and we will be able 
to translate some of the results there to the problem at 
hand. We find, as in[2^ a rich phase diagram, with dif- 
ferent fixed points. The stability of these fixed points 



depends on the nature of the disorder. Finally, we dis- 
cuss the changes introduced in this picture by the possible 
existence of localized states near the Fermi edge induced 
by strong distortions of the lattice. 

The model: Coulomb interaction and disorder. We de- 
scribe the electronic states within each graphene plane 
by two two-component spinors associated to the two in- 
equivalent Fermi points in the Brillouin Zone. They are 
combined to a four component (Dirac) spinor. These 
spinors obey the massless Dirac equation. The Hamilto- 
nian of the free system is: 



(1) 



where "if = \I'^7o with the 4x4 matrix 70 = (T3 (8> (T3. 
We further have 7 = (71,72) = (— «C2,«ci) C?> us. The 
denote the usual Pauli matrices such that {7^1,71/} = 
2(7^,,yl4x4, gp,^i, denoting the Minkowski tensor where 
.90,0 — 1; 9i,i — ~1 with i — 1,2, and zero otherwise. 

The long range Coulomb interaction in terms of the 
Dirac spinors reads 

Hee^^-f I d?xd^x'^{x)^^^{x)—^;^{x')^^^{x') 

An J \x — x'\ 

(2) 

where g = /vp is the dimensionless coupling constant. 

In order to describe disorder effects, the Dirac spinors 
are coupled to a gauge field A{x), 



H, 
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d^x'^{x)V^{x)A{x) 



(3) 



where vy characterizes the strength and the 4x4 matrix 
r the type of the vertex. In general, A{x) is a quenched, 
Gaussian variable with the dimensionless variance A, i.e.. 



{A{x)) = , {A{x)A{x')) = A6^{x-x') 



(4) 



We will discuss five different types of disorder which are 
associated to the five mutually anticommuting 4x4 ma- 
trices plus the unity matrix: i) For a random chemical 
potential, the 4x4 matrix F is given by F = 70. The 
long range components of this type of disorder do not 
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induce transitions between the two inequivalent Fermi 
points. This type of disorder will yield an unstable fixed 
line, ii) A random gauge potential involves the 4x4 
matrices T — iji and T = ij2- This type of disorder 
will yield a stable fixed line which is linear in the {g, A)- 
plane. iii) (a) A fluctuating mass term is described by 
r = 14x4- (b) Topological disorder is given by F = ^75 
with 75 = l2x2 ® o'2- This type of disorder is associated 
to the existence of pentagons and heptagons, or, more 
generally, to local distortions of the lattice axesp^ Il9l|. 
This term can thus be represented by a gauge potential 
which induces transitions between the two Fermi points, 
(c) To complete the discussion, we also mention F = ^75 
where 75 — l2x2 ® cti. This vertex type can be related 
to an imaginary mass that couples the two inequivalent 
Fermi points. All these types of disorder will yield a sta- 
ble fixed line which is cubic in the (g, A)-plane. 

Renormalization of the effective couplings. To discuss 
the renormalizability of the theory, we will first treat 
the disorder gauge-field as an external potential and do 
not consider the average over different realizations of this 
field. The free, massless Dirac propagator is given bv[ll|: 



joLu — VFJ ■ p + iO 

Within one loop and without averaging over the dis- 
order potential, only two diagrams need to be consid- 
ered, i.e., the self-energy of the fermion propagator due 
to electron-electron interaction and the vertex correction 
of the "external" gauge field that couples to the Dirac 
bilinear ^F^*. Both diagrams are shown on the left hand 
side of Fig. p. 

The Fermi velocity is renormalized by the top left dia- 
gram of Fig. ^ aswF = ZypVp with Z^p — 1 — o/(167re), 
using dimensional regularization with e — > 0|ll|. Notice 
that the vertex of the Coulomb interaction is not renor- 
malized within one-loop order. 

Whether or not vr needs to be renormalized depends 
on the type of the disorder: i) For a random chemical 
potential with F = 70, the right diagram at the top 
of Fig. pP is non-divergent. This statement is equiva- 
lent to the one that there is no vertex correction of the 
Coulomb interaction. We thus set = vi with the flow- 
invariant velocity vi. ii) For a random gauge potential 
(F — 171, «72), the vertex is renormalized by the same 
factor as the Fermi velocity, i.e., Zr = 1 — (?/(167re). 
This fact is related to the conservation of the current. 
We can thus set vr = vp- iii) For a random mass term 
F = l4x4, topological disorder F — 275, and F = ^75, the 
vertex strength vr is renormalized by = 1 — g/(87re) 
which follows from i) and ii). Within one loop- level, we 
can thus set vr = Vp/v^ with the flow- invariant velocity 

Notice that - irrespective of the type of disorder - our 
model depends on a single renormalized parameter, vp- 




FIG. 1: Top: Self energy corrections. Dashed line, Coulomb 
interaction. Dotted line, ph external disorder. Bottom: 
Renormalization of the vertex for the external disorder. 
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FIG. 2: One-loop phase diagram for two-dimensional mass- 
less Dirac spinors including long-ranged electron-electron in- 
teraction g and disorder A. Top: Random chemical potential 
(r = 70). Center: Random gauge potential (F = 171,172). 
Bottom: Random mass term (F = 14x4), topological disor- 
der (r = 175), and F = 175. 

Averaging over disorder. As it was discussed in Ref. 
ITlL there is no wave function renormalization Z^ to lead- 
ing order in g. Including self energy corrections due to 
averaging over an ensemble of various realizations of the 
gauge field A{x), the wave function gets renormalized 19]. 
The self energy due to disorder is shown at the top right 
of Fig. and reads 

I]r(^) = 1 ^FGo(c.,K)F . (6) 

The wave function renormalization is independent of the 
vertex type and yields Z^, = 1 — Ai;p/(327rwp£), where 
again we use dimensional regularization with £ ^ 0. 

The wave function renormalization Z$ also changes the 
renormalization factor of the Fermi velocity as we keep 
Gq ^ — S invariant with 

Go^zi^^P) = -Z^{-ioU} - ZypVF'i -P) ■ (7) 
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We also need to discuss the vertex corrections due to 
averaging over the disorder. The diagram which renor- 
mahzes the vertex of the external gauge field is shown at 
the bottom of Fig. Also the vertex of the Coulomb 
interaction is being renormalized (not shown). The ver- 
tex correction to T due to averaging over the disorder 
type r' is given by 



16 



(8) 



This expression depends on the given combination of ver- 
tices, but does not result in a renormalization of the elec- 
tric charge (F = 70). In the case, where there is only 
one type of disorder, the vertex correction exactly com- 
pensates the effect of on Zy^ such that the relation 
between vr and vp remains valid, i.e., i) vr = wi, ii) 
Vr = vp, and iii) vr — Vp/v^ with wi, constant. 

Phase diagrams. Disorder thus only changes the flow 
of the Fermi velocity due to wave function renormal- 
ization, i.e., Zyp —>■ Zyp/Zq,. From the beta- function 
/3yp = Ad\ZypVF, we obtain the following flow equation 
for the effective Fermi velocity v'^f'^ (i' = In A/Aq ^ 1/e): 



d V 



d£ Vp 



1 

16^ 




(9) 



We can now discuss the phase diagram for the vari- 
ous types of disorder: i) For a random chemical potential 
(T = Jo), Vr — vi remains constant under renormaliza- 
tion group transformation. There is thus an unstable 
fixed line at Vp — vfA/{2e^). In the (5, A) -plane, the 
strong-coupling and the weak-coupling phases are sep- 
arated by a hyperbola, with the critical electron inter- 
action g* = e^/vp = 2e'^/{v\/S). ii) A random gauge 
potential involves the vertices F = i"fi,i"f2- The vertex 
strength renormalizes as vr — vp. There is thus an at- 
tractive Luttingcr-like fixed point for each disorder corre- 
lation strength A given by Vp = 2e^/A or g* — A/2, iii) 
For a random mass term F = l4x4, topological disorder 
r = «75, and F = 175, we have vr — Vp/v^. There is thus 
again an attractive Luttinger-like fixed point for each dis- 
order correlation strength A given by Vp = \j2v\e'^ j A. 

or 5* = </Ae4/(2«i)- _ 

To make connection to previous work [2^, we define 

A = ^{v!^ ^ jv'^p^)'^ . This yields the linear fixed line 
g* ~ A with i) A cx 5^, ii) A cx const., and iii) A cx g^^ . 
For one disorder type, our results thus agree with the 
ones by Ye. 

Localized states. The tight binding model defined 
by the tt orbitals at the lattice sites can have edge 
states when the sites at the edge belong all to the same 
sublatticeli^, 0, H^- These states lie at zero energy, 
which, for neutral graphene planes, correspond to the 
Fermi energy. The states at zero energy are localized in 
one of the two interpenetrating sublattices which can be 




f(z)=0 



FIG. 3: Top: example of a crack in a graphene plane. The 
atoms at the upper edge and those at the lower edge belong 
to different sublattices. Bottom: approximate cut in the com- 
plex plane which can be used to represent this crack at long 
distances. See text for details. 



defined in the honeycomb structure. In terms of each of 
the two component spinors, in the continuum approxi- 
mation used here, a possible state localized in one edge 
can be written as ^'(r) = A{e^'^^ , 0), where z = x + iy, we 
assume that the edge lies at y = 0, and that the graphene 
fills the upper half plane, y > 0. 

In a strongly disordered sample, large defects made up 
of many vacancies can exist. These defects give rise to 
localized states, when the termination at the edges is lo- 
cally similar to the surfaces discussed above. Note that, if 
the bonds at the edges are saturated by bonding to other 
elements, like hydrogen, the states at these sites are re- 
moved from the Fermi energy, but a similar boundary 
problem arises for the remaining tt orbitals. A particular 
simple example is given by the crack shown in Fig.[3|. 
The boundary conditions are such that ^'(2) = [/(2:),0] 
where f{z) — at the surface of the crack, because it has 
edges of the two sublattice types. Possible edge states 



are: \E'„ 



Re 



A similar solution is 



obtained by exchanging the upper by the lower spinor 
component, and replacing z ^ z. Because of the dis- 
creteness of the lattice, the values of n should be smaller 
than the number of lattice units spanned by the crack. 

These states are half filled in a neutral graphene plane. 
In the absence of electron interactions, this leads to a 
large degeneracy in the ground state. A finite local re- 
pulsion will tend to induce a ferromagnetic alignment of 
the electrons occupying these states, as in similar cases 
with degenerate bands |2a|. Hence, we can assume that 
the presence of these states leads to magnetic moments 
localized near the defects. 

We now have to analyze the influence of these magnetic 
moments in conduction band described in the previous 
sections. The hopping between the states involved in the 
formation of these moments and the delocalized states in 
the conduction band vanishes when the localized states 
lie at zero energy. Hence, an antifcrromagnetic Kondo 
coupling will not be induced. The localized and conduc- 
tion band states, on the other hand, coexist at the same 
lattice sites. The presence of a finite local repulsion, U, 
will lead to a ferromagnetic coupling. The local coupling, 
at site i, between the localized states and the conduction 
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band is proportional to UJ^j Pi,j' where pij is charge of 
state j at site i. The conduction electrons will mediate 
an RKKY interaction between the localized moments: 

jRKKvir) ^U' f d'ke^*^'^x(k) ^ C/'^^ (10) 
J VF\r\-^ 

Where the static susceptibility is x(k) oc |k|[23, and a is 
the lattice constant. It is interesting to note that, due to 
the absence of a finite Fermi surface, the RKKY interac- 
tion in ea. H10() does not have oscillations. Hence, there 
are no competing ferro- and antiferromagnetic couplings, 
and the magnetic moments will tend to be ferromagnet- 
ically aligned. Asuming U ^ hvp/a ^ 1 — 3eV, this 
argument leads to a Curie temperature Tc ^ U{a^/d'^), 
where d is the distance between defects. 

Conclusions. In this work, we estimated the effect of 
disorder in graphene sheets in the presence of long-ranged 
electron-electron interaction. The disorder term can be 
characterized by five different vertices which lead - within 
a one-loop level - to three distinct phase diagrams. A ran- 
dom chemical potential divides the phase diagram into a 
strong-coupling and weak-coupling regime by an unstable 
fixed line. A random gauge potential as well as a random 
mass term and topological disorder alter the electronic 
structure of two-dimensional graphene sheets as it drives 



the system towards a stable, Luttinger-like fixed point. 
As extensively analyzed previously (see, for instance jljj) 
the dimensionless Coulomb coupling, for graphene planes 
is g = e^/^egvp) 1. Randomness in the chemical po- 
tential, due to local defects, like impurities or vacancies, 
lead to A {V/fiVF)'^c, where V is the strength of the 
local potential and c is the concentration (A ~ ca^ for va- 
cancies, where a is the lattice constant). Lattice defects, 
such as inclinations or disclinations, induce an effective 
gauge potential, whose strength, for dislocations, can be 
expressed as A ~ cb^ 19], where b is the Burgers vector of 
the dislocation. The actual strength of the disorder, how- 
ever, can be much larger for experiments performed on 
irradiated samples p3|. Hence, it is likely that graphene 
sheets are in the intermediate range coupling, and we 
can expect a qualitative behavior like that reported here. 
Note that the presence of disorder at intermediate scales 
will suppress the chiral symmetry breaking transition ex- 
pected for pure graphene 13, 14] . 
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